Antiferromagnetic Spin Fluctuations and the Pseudogap in the Paramagnetic Phases 
of Quasi-Two-Dimensional Organic Superconductors 
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The metallic states of a broad range of strongly correlated electron materials exhibit the subtle 
interplay between antiferromagnetic spin fluctuations, a pseudogap in the excitation spectra, and 
non-Fermi liquid properties. In order to understand these issues better in the k-(ET)2X family 
of organic charge transfer salts we give a quantitative analysis of the published results of nuclear 
magnetic resonance (NMR) experiments. The temperature dependence of the nuclear spin relaxation 
rate l/Ti, the Knight shift Ks, and the Korringa ratio KL, are compared to the predictions of 
the phenomenological spin fluctuation model of Moriya, and Millis, Monien and Pines (M-MMP), 
that has been used extensively to quantify antiferromagnetic spin fluctuations in the cuprates. For 
temperatures above Tnmr — 50 K, the model gives a good quantitative description of the data for the 
paramagnetic metallic phase of several k-(ET)2X materials, with an antiferromagnetic correlation 
length which increases with decreasing temperature; growing to several lattice constants by Tnmr- It 
is shown that the fact that the dimensionless Korringa ratio is much larger than unity is inconsistent 
with a broad class of theoretical models (such as dynamical mean-field theory) which neglect spatial 
correlations and/or vertex corrections. For materials close to the Mott insulating phase the nuclear 
spin relaxation rate, the Knight shift and the Korringa ratio all decrease significantly with decreasing 
temperature below Tnmr- This cannot be described by the M-MMP model and the most natural 
explanation is that a pseudogap opens up in the density of states below Tnmr, as in, for example, the 
underdoped cuprate superconductors. An analysis of the Mott insulating phase of k-(ET)2Cu2(CN)3 
is somewhat more ambiguous; nevertheless it suggests that the antiferromagnetic correlation length 
is less than a lattice constant, consistent with a large frustration of antiferromagnetic interactions 
as is believed to occur in this material. We show that the NMR measurements reported for k- 
(ET)2Cu2(CN)3 are qualitatively inconsistent with this material having a ground state with long 
range magnetic order. A pseudogap in the metallic state of organic superconductors is an important 
prediction of the resonating valence bond theory of superconductivity. Understanding the nature, 
origin, and momentum dependence of the pseudogap and its relationship to superconductivity are 
important outstanding problems. We propose specific new experiments on organic superconductors 
to elucidate these issues. Specifically, measurements should be performed to see if high magnetic 
fields or high pressures can be used to close the pseudogap. 
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I. INTRODUCTION 



In the past twenty years a diverse range of new strongly 
correlated electron materials with exotic electronic and 
magnetic properties have been synthesized. Examples in- 
clude high-temperature cuprate superconductors,^ man- 
ganites with colossal magnetoresistance^ cerium ox- 
ide catalysts,^ sodium cobaltates,™ ruthenates^"^ heavy 
fermion materials,^ and superconducting organic charge 
transfer saltsi^ Many of these materials exhibit a sub- 
tle competition between diverse phases: paramagnetic, 
superconducting, insulating, and the different types of 
order associated with charge, spin, orbital, and lattice 
degrees of freedom. These different phases can be ex- 
plored by varying experimental control parameters such 
as temperature, pressure, magnetic field, and chemical 
composition. Although chemically and structurally di- 
verse the properties of these materials are determined 
by some common features; such as, strong interactions 
between the electrons, reduced dimensionality associated 
with a layered crystal structure, large quantum fluctu- 
ations, and competing interactions. Many of these ma- 



terials are characterized by large antiferromagnetic spin 
fluctuations. Nuclear magnetic resonance spectroscopy 
has proven to be a powerful probe of local spin dynamics 
in many strongly correlated electron materialsi^ii^iiiii^ 
Longer range and faster spin dynamics have been studied 
with inelastic neutron scattering. One poorly understood 
property of the paramagnetic phases of many of these 
materials is the pseudogap present in large regions of the 
phase diagram. Although the pseudogap has received 
the most attention in the cuprateSj^ it is also present in 
quasi one-dimensional charge-density wave compoundsji^ 
manganitesr^ heavy fermion materials ji^, in simple met- 
als with no signs of superconductivity,!^ and quite pos- 
sibly in organic charge transfer materials. The focus of 
this paper is on understanding what information about 
antiferromagnetic spin fluctuations and the pseudogap 
can be extracted from NMR experiments on the organic 
charge transfer salts. 

The systems which are the subject of the current 
study are the organic charge transfer salts based on 
electron donor molecules BEDT-TTF (ET), in particu- 
lar the family k-(ET)2X (where k indicates a particular 
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polymorph''^^). Remarkably similar physics occurs in the 
other dimerised polymorphs, such as the /3, /?', A, and k 
phases.— These materials display a wide variety of uncon- 
ventional behaviours^ including: antiferromagnetic and 
spin liquid insulating states, unconventional supercon- 
ductivity, and the paramagnetic metallic phases which we 
focus on in this paper. They also share highly anisotropic 
crystal and band structures. However, as for various so- 
ciological and historical reasons, the k salts have been 
far more extensively studied, and because we intend, in 
this paper, to make detailed comparisons with experi- 
mental data, we limit our study to k phase salts. This of 
course begs the question do similar phenomena to those 
described below occur in the f3, (3' , oi X salts? We would 
suggest that the answer is probably yes but this remains 
an inviting experimental question. 

As well as their interesting phenomenology the organic 
charge transfer salts are important model systems and a 
deeper understanding of these materials may help address 
a number of important fundamental questions concerning 
strongly correlated systems. The (non-interacting) band 
structure of the k phase salts is well approximated by the 
half-filled tight binding model on the anisotropic lattice.'^ 
This model has two parameters t, the hopping integral 
between nearest neighbor (ET)2 dimers, and t' , the hop- 
ping integral between next nearest neighbors across one 
diagonal only. In order to describe strongly correlated 
phases of these materials we must also include the effects 
of the Coulomb interaction between electrons. The sim- 
plest model which can include these strong correlations 
in the Hubbard model contains one additional parameter 
over the tight binding model: U, the Coulomb repulsion 
between two conduction electrons on the same dimcr. In 
the Hubbard model picture the different k-{F,T)2X salts 
and different pressures correspond to different values of 
t' /t and U/t, but all of the k phase salts are half filled. 
Varying U/t allows us to tune the proximity to the Mott 
transition - understanding the Mott transition and its as- 
sociated phenomena remains one most important prob- 
lems in theoretical physics22iS^ and the organics have pro- 
vided a new window on this problem4i2ii2^i2£i^ Varying 
t' /t allows one to tune the degree of frustration in the sys- 
tem. Understanding the effect of frustration in strongly 
correlated systems is of general importance. "'■^^'^^ For 
example there are strong analogies between the organ- 
ics and Naa;Co02,'^° and much recent interest has been 
sparked by the observation of possible spin liquid states 
in organic charge transfer salts i^i^i^^ 

The paramagnetic metallic phases of k-(ET)2^ are 
very different from a conventional metallic phase. Many 
features of the paramagnetic metallic phases agree well 
with the predictions of dynamical mean field theory 
(DMFT) which describes crossover from 'bad metal' at 
high temperatures to a Fermi liquid as the tempera- 
ture is lowered j^^i^i^i^ This crossover from incoher- 
ent to coherent intralayer-i^^ transport has been ob- 
served in a number of experiments such as resistivity^^^ 
thermopower ;^^!^^ and ultrasonic attenuationi^i^ The 



existence of coherent quasiparticles is also apparent from 
the observed magnetic quantum oscillations at low tem- 
peratures in K-(ET)2^j^'^2iiL However, nuclear mag- 
netic resonance experiments (reproduced in Figs. [T]and 
[2]) on the paramagnetic metallic phases on k-(ET)2^ do 
not find the well known properties of a Fermi liquid. 
The nuclear spin relaxation rate per unit temperature, 
1/TiT, is larger than the Korringa form predicted from 
Fermi liquid theory. As the temperature is lowered 
l/TiT reaches a maximum; we label this temperature 
Tnmr (the exact value of Tnmr varies with the anion 
X, but typically, Tnmr ~ 50 K, see Fig. O. l/TiT de- 
creases rapidly as the temperature is lowered below Tnmr 
[see Fig [i] . 12,18,19 r^-^^ Knight shift also drops rapidly 
around Tnmr [see Fig|4],^^ This is clearly in contrast the 
Korringa-like behavior one would expect for a Fermi liq- 
uid in which l/TiT and are constant for T <C Tp, 
the Fermi temperature. Similar non-Fermi liquid tem- 
perature dependences of 1/TiT and are observed in 
the cuprates42ii^ For the cuprates, it has been argued 
that the large enhancement of l/TiT is associated with 
the growth of antiferromagnetic spin fluctuation within 
the Cu02 planes as the temperature is lowered4^i^ The 
large decrease observed in l/TiT and Kg is suggestive of 
a depletion of the density of states (DOS) at the Fermi 
level which might be expected if a pseudogap opens at 
Tnmr- 

A qualitative description of spin fluctuations in the 
paramagnetic metallic phase of k-(EiT)2X has not been 
performed previously. The importance of spin fluctua- 
tions in K-(ET)2X below Tc has been pointed out by 
several groups . ^''^^I'^^'^^i'*^'^" Since the nature of the para- 
magnetic metallic and superconducting phases are inti- 
mately intertwined (in most theories, including BCS, su- 
perconductivity arises from an instability of the metal- 
lic phase), it is important to understand whether the 
spin fluctuations may extend beyond the superconduct- 
ing region into the paramagnetic metallic phase and if 
so how strong they are. Another unresolved puzzle is 
whether there is a pseudogap in the paramagnetic metal- 
lic phase of k-(ET)2X which suppresses the density of 
states (DOS) at the Fermi level. A pseudogap has been 
suggested on the basis of the NMR data and heat ca- 
pacity measurement.^^ If there is a pseudogap then im- 
portant questions to answer include: (i) how similar is 
the pseudogap in k-(ET)2X to the pseudogaps in the 
cuprates, manganites and heavy fermions? (ii) is the 
pseudogap in k-(ET)2X related to superconductivity? 
and (iii) if so how? 

A number of scenarios in which a pseudogap may arise 
have been proposed. One possible origin of a pseudogap, 
which a number of authors^'^'^i^i^ have argued may 
be relevant to the organic charge transfer salts, is the 
resonating valence bond (RVB) picture (for a review see 
Refs. [H and Issl) . In this picture the electron spins form 
a linear superposition of spin singlet pairs. The singlet 
formation can naturally explain the appearance of a gap: 
a non-zero amount of energy is required to break the sin- 
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FIG. 1: [Color online] Comparison of the measured nuclear spin relaxation rate per unit temperature, 1/T\T, with the predic- 
tions of the spin fluctuation model for various organic charge transfer salts. Panel (a) shows data for K-(ET)2Cu[N(CN)2]Br mea- 
sured by Mayaffre et al}^ , panel (b) shows data for «:-(ET)2Cu[N(CN)2]Br measured by De Soto et ali^, panel (c) shows data 
for K(d8)-(ET)2Cu[N(CN)2]Br measured by Miyagawa et ali^, and panel (d) shows data for a fi;-(ET)2Cu(NCS)2 powder 
sample measured by Kawamoto et al.S^ The 1/TiT data are weakly temperature dependent at high temperatures, have a 
maximum at Tnmr ~ 50 K, and drop abruptly below Tnmr, contrary to what one would expect for a Fermi liquid in which 
1/T\T is constant. The remarkable similarities of these data results from the quantitative and qualitative similarity of the 
antiferromagnetic spin fluctuations in the paramagnetic metallic phases of these materials. The parameters that produce the 
best flts (solid lines) to Eq. (|13p are tabulated in Table |I] The spin fluctuation model gives a good flt to the experimental 
data between Tnmr ~ 50 and room temperature which suggests strong spin fluctuations in the paramagnetic metallic states of 
K-(ET)2Cu[N(CN)2]Br, fi:(d8)-(ET)2Cu[N(CN)2]Br, and k-(ET)2Cu(NCS)2. However, below Tnmr the spin fluctuation model 
does not describe the data well. In section |V] we argue that this is because a pseudogap opens at Tnmr. In each figure the 
solid line indicates the large T approximation of the spin fluctuation model [Eq. (|15b|) ]. To check that this approximation 
is reasonable we also plot the full spin fluctuation model [Eq. (|13|) ] as a dashed line in panel (b). The full and dashed lines 
cannot be distinguished until well below Tnmr and so we concluded that the high T approximation is excellent in the relevant 
regime. Note that the analysis on 1/T\T cannot differentiate between antiferromagnetic and ferromagnetic spin fluctuations 
(see section fll A 2|) . but the Korringa ratio strongly differentiates between these two case and indicates that the fluctuations 
are antiferromagnetic (see Fig. [2]). The nomenclature K-Br, d8-Br, and k-NCS is used as shorthand for K-(ET)2Cu[N(CN)2]Br, 
K(d8)-(ET)2Cu[N(CN)2]Br, and k-(ET)2Cu(NCS)2 respectively in the figure keys. 



glet pairs. For weakly frustrated lattices, such as the 
anisotropic triangular lattice in the parameter range ap- 
propriate for K-(ET)2Cu[N(CN)2]Br, k-(ET)2Cu(NCS)2, 
and K-(ET)2Cu[N(CN)2]Cl, the RVB theory predicts 'd- 
wave' superconductivity^ii^i^i^'^iS This is the symme- 
try most consistent with a range of experiments on these 
K-(ET)2X saltsJ^i^i^i^ However, as the frustration is 
increased changes in the nature of the spin fluctuations 
drive changes in the symmetry of the superconducting 
state<^i^ For example, for the isotropic triangular lat- 
tice t = t' [t ^ t' is thought to be appropriate for n- 
(ET)2Cu2(CN)3] RVB theory predicts that the super- 



conducting order parameter has 'd -I- id' symmetry , ^'^'^-'^ 

RVB theory also predicts that a pseudogap with the 
same symmetry as the superconducting state exists in 
the paramagnetic phase at temperatures above the su- 
perconducting critical temperature, Tc- This pseudogap 
results from the formation of short range singlets above 
T(,\ only at Tc do these singlets acquire off-diagonal long- 
range order. There are two energy scales (A and A in 
the notation of, e.g., Refs. \^ and IsgI) in the RVB the- 
ory. In the simplest reading of the theory the ratio 
A/A « r*/Tc where T* is the temperature at which the 
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FIG. 2: [Color online] Comparison of the Korringa ratio 
K. (X l/TiTKl of K-(ET)2Cu[N(CN)2]Br measured by De 
Soto et al^ with the prediction of the antiferromagnetic spin 
fluctuation model. The best fit to Eq. (|16|l is indicated by 
the solid line. The Korringa ratio is larger than 1 which indi- 
cates that the spin fluctuations are antiferromagnetic (/C < 1 
for ferromagnetic fluctuations, see Section mXlJ. The an- 
tiferromagnetic correlation length is found to be 3.5 ± 2.5 
lattice spacings at 50 K. Below 50 K the Korringa ratio is 
suppressed, and in section |V] we argue that this is because 
a pseudogap opens at Tnmr. In the key to this flgure k- 
(ET)2Cu[N(CN)2]Br is abbreviated as K-Br. 



pseudogap opens. For the appropriate model Hamilto- 
nian for the layered hi-{W£)2X salts A/A ~ 5 near the 
Mott transition^ which is remarkably similar to the ratio 
TNMR/Te in K-(ET)2Cu[N(CN)2]Br. 

We use the phenomenological antiferromagnetic spin 
fluctuation model which was first introduced by Moriya^ 
and then applied by Millis, Monien and Pines (MMP)i^ 
to cuprates, to examine the role of spin fluctuations in the 
paramagnetic metallic phase of k-(ET)2X. We investi- 
gate whether the anomalous temperature dependences of 
NMR data can be explained without invoking a pseudo- 
gap in the DOS. We fit the spin fluctuation model to the 
nuclear spin relaxation rate per unit temperature 1/TiT, 
Knight shift Ks, and Korringa ratio K. data and find that 
the large enhancements measured in 1/TiT and K. above 
Tnmr are the result of large antiferromagnetic spin fluc- 
tuations [see Figs. [T]and[2]. The correlation of the anti- 
ferromagnetic spin fluctuations increases as temperature 
decreases and the relevant correlation length is found to 
be 3.5 ± 2.5 lattice spacings at T = 50 K. The model 
produces reasonable agreement with experimental data 
down to T - 50 K. Below 50 K, l/TiT is suppressed 
but never saturates to a constant expected for a Fermi 
liquid while the spin fluctuation model produces a mono- 
tonically increasing l/TiT with decreasing temperature. 
We argue that this discrepancy between theory and ex- 
periment is due to the appearance of a pseudogap, not 
captured by the spin fluctuation model, which suppresses 



the DOS at the Fermi level. 



Our results suggest the paramagnetic metallic phase 
of K-(ET)2^ is richer than a renormalized Fermi liquid 
as has been previously thought to describe the low tem- 
perature metallic state. An exotic regime similar to the 
pseudogap in the cuprates appears to be realized in the 
paramagnetic metallic phase of k-(ET)2X . Thus we be- 
lieve the appropriate phase diagram of k-(ET)2^ looks 
like the one sketched in Fig. [H However, the pseudogap 
in K-(ET)2X is rather peculiar. On one hand, it shows a 
coherent intralayer transport (apparent from the re- 
sistivity and the observed quantum oscillations). On the 
other hand it also shows a loss of DOS apparent from 
1/TiT and Kg. Therefore it is important to emphasize 
that the pseudogap phase proposed for k-{ET)2X is dif- 
ferent from the pseudogap phase realized in the cuprates. 
Understanding these differences may well provide impor- 
tant insight into the physics of the pseudogaps in both 
classes of material. We will discuss this matter further 
in Section V. 



We have also applied the spin fluctuation formalism 
to the Mott insulating phase k-(ET)2Cu2(CN)3 (which 
may have a spin liquid ground state). '^^ While a rea- 
sonable agreement between the calculated and the ex- 
perimental data on 1/TiT has been obtained, the re- 
sult for the Knight shift does not agree with the data. 
We believe that this discrepancy reflects the failure of 
the assumption that the peak in the dynamic suscepti- 
bility dominates even the long wavelength physics im- 
plicit in the spin-fluctuation model. The failure of this 
assumption is consistent with the strong frustration in 
k-(ET)2Cu2(CN)3 and the observed spin-liquid behavior 
of this material.-S»31 We show that there are qualitative 
as well as quantitative difference between the spin fluc- 
tuations in k-(ET)2Cu2(CN)3 and the spin fluctuations 
in the other K-phase salts discussed in this paper. 



The structure of the paper is as follows. In Section II 
we introduce the temperature dependence of the nuclear 
spin relaxation rate, spin echo decay rate. Knight shift, 
and Korringa ratio and describe how they probe spin 
fluctuations by probing the dynamical susceptibility. We 
calculate these properties in a number of approximations 
and contrast the results. In Section III we demonstrate 
that the spin fluctuation model provides reasonable flts to 
the existing experimental results for k-(ET)2X and dis- 
cuss its limitations when applied to those materials. Sec- 
tion IV deals with the application of the spin fluctuation 
model to the spin hquid compound k-(ET)2Cu2(CN)3. 
We discuss the nature and extent of the pseudogap in 
Section V. Finally, in Section VI we suggest new experi- 
ments and give our conclusions. 
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FIG. 3: [Color online] The schematic phase diagram for k-(ET)2X as a function of temperature and pressure. Thin solid lines 
represent second order phase transitions, the thick solid line is the first order transition line which ends at a critical point shown 
as a filled circle, and dashed lines indicate crossovers. The pseudogap phase is much more complicated than a renormalized 
Fermi liquid that has been previously thought to characterize the paramagnetic metallic phase at low temperatures: it shows 
a coherent transport character with long lived quasiparticles, marked by resistivity behaviop^^ and magnetic quantum 
oscillations,^^ but at the same time it exhibits a loss of density of states which is clearly seen in the NMR data. There are 
not sufficient data at this moment to determine where the pseudogap phase boundary ends so this uncertainty is represented 
by the shaded area with the question mark. More detailed experimental and theoretical studies in the vicinity of this point 
will give important insight into how the pseudogap is related to superconductivity. The possibility of a quantum critical 
point somewhere in the vicinity of the point where the superconducting critical temperature goes to zero may have important 
consequences for the observation that the materials with the lowest superconducting critical temperatures have extremely small 
superfluid stiffnesses and are very different from BCS superconductors.— i^^*^ 



II. THE SPIN LATTICE RELAXATION RATE, 
SPIN ECHO DECAY RATE, AND KNIGHT 
SHIFT 



The temperature dependence of the nuclear spin lattice 
relaxation rate l/TiT, spin echo decay rate I/T2, Knight 
shift Ks, Korringa ratio /C, and the real and imaginary 
part of the dynamic susceptibility, x'i^i ^) ^-^d x"(<3i 
defined by 



X(q,w)= / ^^te*('^■"-'^*)x(I■,^)■ 



(l) 



are discussed in this section. The general expressions for 
these quantities are^ 



2X"(q,^) 



^ - ||El^(q)lV(q,0)^ 



and 



l^(o)lx'(o,o) ^ 



Airks \"fN J TiTKl 



(2a) 
(2b) 
(2c) 

(2d) 



where A{c\) is the hyperfine coupling between the nuclear 
and electron spins, 7Ar (7e) is the nuclear (electronic) gy- 
romagnetic ratio, and /a is the relative abundance of 
the nuclear spin. For simplicity we will often consider a 
momentum independent hyperfine coupling |^| in what 
follows. Note that Eqs. [5] show that this is an uncon- 
trolled approximation for both Ti and T2, but that it is 
not an approximation at all for Kg- This is because Ks 
only probes the long wavelength physics and hence only 
depends on A(0), the hyperfine coupling at q = 0. 



The calculation of the quantities in Eqs. ^ boils down 
to determining the appropriate form of the dynamic sus- 
ceptibility. In the following sections we begin by dis- 
cussing the role of vertex corrections in determining the 
properties measured by NMR (section ??), before moving 
onto the a variety of approximations for calculating the 
dynamic susceptibility. They are antiferromagnetic and 
ferromagnetic spin fluctuations fscction lll A[) . dynamical 
mean field theory fsection FlIBp . and the approach to 
the quantum critical region of frustrated two-dimensional 
antiferromagnets fsection lll C\ . However, in this work we 
will predominantly use the spin fluctuation model to an- 
alyze the NMR data. The other models are presented for 
comparison. 
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A. The Spin Fluctuation Model 



The dynamic susceptibility in this model is written 

„44.45 



ha! 



,45 



x(q, ^) = XLw(t^) + XAF(q, w), 



(3) 



where Xhwi^) is the dynamic susceptibility in the long 
wavelength regime and xaf(<1i '-^) is a contribution to the 
dynamic susceptibility which peaks at some wave vector 
Q. These susceptibilities take the form 



XAF(q,t^) 



Xo{T) 



1 - iuj/T(T) 



xq{t) 



i + e(T)'|q-Qp- W^sF(r) 



(4) 



where Xo{T) [xq(2^)] is the static spin susceptibility at 
q = [Q], r(r) [wsfC?")] is the characteristic spin fluc- 
tuation energy which represents damping in the system 
near q = [Q], and ^(T) is the temperature dependent 
correlation length. Hence, the real and imaginary parts 
of the dynamic susceptibility can then be written as 



x'(q,o) 

X"(q,^) 



Xo(T) 
^Xa{T) 

1 



1 



Xq(T) 



1 



Xo(T) (l + e(T)2|q-Q| 



xq{t)t 



xo(rVsF(r) (i + e(r)2|q-Q|2)2 



(5) 



Note that the above form ol xlw{^) is the appropriate 
form for a Fermi liquid. Therefore if the system under 
discussion is not a Fermi liquid then the validity this ex- 
pression for Xlw{^) cannot be guaranteed. For example, 
the marginal Fermi liquid theory predicts a different fre- 
quency dependence 1^ If the dynamic susceptibility is suf- 
ficiently peaked then 1 /Ti will not be strongly dependent 
on the long wavelength physics [because 1/Ti measures 
the susceptibility over the entire Brillouin zone, c.f., Eq. 
(l2a| . and therefore will be dominated by the physics at 
q = Q]. On the other hand the Knight shift is a mea- 
sure of the long wavelength properties [c.f., Eq. (Pc)) ] 
and therefore may be sensitive to the details of xlw(^)- 
Here we will explicitly assume, following the assumption 
made by MMP"*^, that the uniform susceptibility xo and 
the spin fluctuation energy near g = is temperature 
independent. One justification for this approximation 
in organics is that the Knight shift is not strongly tem- 
perature dependent [see Section HI D]. This approxima- 
tion will break down in the systems where the uniform 
susceptibility is strongly temperature dependent such as 
YBa2Cu306.63^ and Lai.8Sro.i5Cu04..^ 

In the critical region ^(T) ^ a, where £,{T) is the cor- 
relation length, and a is the lattice constant, and one 



2-r, 



xq{t) - 



Xo 



(6) 



where rj is the critical exponent which governs the power- 
law decay of the spin correlation function, z is the dy- 
namical critical exponent, and is some temperature 
independent length scale. The simplest possible assump- 
tions are relaxation dynamics for the spin fluctuations 
(which are characterized by z = 2) and mean fleld scal- 
ing of the spin correlations {rj ~ 0). With these approx- 
imations the real and imaginary parts of the dynamic 
susceptibility are given by 



X'(q,0) = 



X"(q,^) = 



Xo 
^Xo 

r 



l + f3 



[l+e(r)2|q-Q| 
K(T)/a]4 



[l+e(r)2|q-Q| 



212 



(7) 



where /3 = (a/^o)^- The temperature independent, di- 
mensionless parameter /3 can also be expressed in terms 
of the original variables appearing in the dynamic sus- 
ceptibility in Eq. ^ as 



XQ{T)r f a 

xo^sF(r) 



(8) 



Written in this form, [3 has a clear interpretation: it rep- 
resents the strength of the spin fluctuations at a finite 
wave vector Q. We will now consider two cases: antifer- 
romagnetic and ferromagnetic spin fluctuations. 



1. Antiferromagnetic Spin Fluctuations 

If we have antiferromagnetic spin fluctuations then the 
dynamic susceptibility x(q, w) peaks at a finite wave vec- 
tor q = Q; for example, for a square lattice Q — (7r,7r). 
The NMR relaxation rate, spin echo decay rate, and 
Knight shift can be calculated straightforwardly from the 
real and imaginary parts of the dynamic susceptibility 
given in Eq. 



1 
1 

^2 



([5]). The results are 
27rfcs|A|'xo 



VP 



1 + 



l+(3- 



4 1 



ln(l 



nil 



1 + /3 



1 + m{T)? 



2rxo 



(9a) 

(9b) 
(9c) 

(9d) 
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where Q is a cutoff from the momentum integration [c.f. 
Eq. m]- For » a: ^ ^T?, l/T^ ^ ^(T), 

and Kg ^ constant which leads to the Korringa ra- 
tio K. ~ {h-y^/2Txo)[Qi{T)]^- In this model the Kor- 
ringa ratio can only be equal to unity if the spin fluc- 
tuations are completely suppressed (/3 = 0). Hence, 
one expects /C > 1 if antiferromagnetic fluctuations are 
dominant fS^'iS This indicates that there are significant 
vertex corrections when there large strong antiferromag- 
netic fluctuations. 



2. Ferromagnetic Spin Fluctuations 

For ferromagnetic spin fluctuations, xil^ ^) is peaked 
at q = 0. The NMR relaxation rate and spin echo decay 
rate are exactly the same as those given in Eqs. ([9a|) 
and (j9b[) because 1/TiT and I/T2 come from summing 
the contributions form all wave vectors in the first Bril- 
louin zone, which makes the location of the peak in q 
space irrelevant. In contrast, the Knight shift will be dif- 
ferent in the ferromagnetic and antiferromagnetic cases 
because Ks only measures the q = part of the dynamic 
susceptibility; it will be enhanced by the ferromagnetic 
fluctuations. Thus, in the ferromagnetic spin fluctuation 
description Kg is given by 



Kg 



\A\xo 



1 + y^(^/a) 



and 



K = 



l + (3 



\aT)laX 



1 + [QC(T)]2 



2rXo [l + ^{^/af 



(10) 



(11) 



For i{T) » a: l/T^T ^ ^{Tf, l/T^ ^ i{T), and 
Kg - C(T)2 which leads to /C ~ {h-il/2Vxo)['n£.{T)/a]-'^. 
We can see that /C < 1 in the presence of ferromagnetic 
fluctuations-Si^ So again vertex corrections are impor- 
tant if the system has strong ferromagnetic fluctuations. 
Recall that, in contrast, in the case of antiferromagnetic 
fluctuations the Korringa ratio is larger than one. Thus 
analyzing the Korringa ratio allows one to straightfor- 
wardly distinguish between antiferromagnetic and ferro- 
magnetic spin fluctuations. 



B. Dynamical Mean Field Theory 



l/TiT - a -h b{To/T) (c.f., Ref. [711) and incoherent 
charge transport. This regime is often refereed to as the 
'bad metal'. Microscopically the bad metal is charac- 
terized by quasi-localized electrons and the absence of 
quasiparticles. Thus DMFT predicts that at high tem- 
peratures l/TiT ^ a + b{To/T), but 1/TiT saturates to a 
constant below Tq. This temperature dependence is sim- 
ilar to that for the single impurity Anderson model.— 
Note that this temperature dependence is qualitatively 
similar to that found for spin fluctuations [c.f., Eq. (|15b[) . 
the high temperature limit of Eq. ([5a|) ] . 

The predictions of DMFT correctly describe the prop- 
erties of a range of transport and thermodynamic exper- 
iments on the organic charge transfer saltsi^ i^"^'^^ This 
suggests that these systems undergo a crossover from a 
bad metal regime for T>Tq to a renormalized Fermi liq- 
uid below Tq. As we shall see in more detail later (see 
Fig. [1]) the nuclear spin relaxation rate is suppressed 
but never saturates below Tnmr; this is not captured by 
DMFT. This suggests that the low-temperature regime 
of K-(ET)2X is more complicated than the renormalized 
Fermi liquid predicted by DMFT which until now, is 
widely believed to be the correct description of the low 
temperature paramagnetic metallic state in the organic 
charge transfer salts. We will discuss the reasons for and 
implications of the failure of DMFT to correctly describe 
NMR experiments on the organic charge salts in section 
lYl 



C. Quantum Critical Region of Frustrated 2D 
Antiferromagnets 

The static uniform and dynamic susceptibilities of 
nearly-critical frustrated 2D antiferromagnets has been 
studied by Chubukov et al^ They considered a long- 
wavelength action with an N component, unit-length, 
complex vector which has a SU(iV)(g)0(2) symmetry and 
performed a 1/iV expansion. This gives susceptibili- 
ties which follow a universal scaling form. As one ap- 
proaches the quantum critical point, the spin stiffness 
will vanish but the ratio between in-plane and out-of- 
plane stiffnesses remains flnite and approaches unity. In 
this regime, to order 1/N, the quantities in Eq. ([2]) are 
given bj*2^ 



1/Ti-T''; I/T2-T''- 



(12) 



DMFT is an approach based on a mapping of the 
Hubbard model onto a self-consistently embedded An- 
derson impurity model.— i^Siill DMFT predicts that the 
metallic phase of the Hubbard model has two regimes 
with a crossover from one to the other at a tempera- 
ture Tq. For T<Tq the system is a renormalized Fermi 
liquid characterized by Korringa-like temperature de- 
pendence of 1/TiT and coherent intralayer transport. 
Above To, the system exhibits anomalous properties with 



III. SPIN FLUCTUATIONS IN k-{ET)2X 

The NMR relaxation rate per unit temperature. 
Knight shift, and Korringa ratio in the antiferromag- 
netic spin fluctuations model are given by Eqs. (Pa)) . 
([5c)) . and ([9d)) . Their temperature dependence comes 
through the antiferromagnetic correlation length. We 
adopt the form of ^(T) from M-MMPii^: ^(r)/^(r^) = 
y/2Tx/{T + Tx). For this form of the correlation length 
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Tx represents a natural temperature scale and ^{T) is 
only weakly temperature dependent for T <^ T^. For 
this choice of ^ (T) /a we have 



1 



1 



1 



1 

K, — /Cn- 



(T/T, + 1)2 + 2^2^(r/r, + 1) 



1 + 27r2C + T/T,j, 

/3C^ 

(T/T, + l)^+27r^C(T/T, + l) 



1 



(13) 



where we have defined 

(i/Tir)o 

and /Co 

to simplify the notation. 



= 2 



n 2 



27rfcB|A|2xo 



2rxo' 



(14) 



A. The Nuclear Spin Relaxation Rate 

We now analyze the temperature dependence of 1 /Ti . 
In the discussion to follow, we will assume that the cor- 
relation length is large compared to unity and to T /Tx^ 
i.e, C = 2(^(r2;)/a)2 > T/T^. By this assumption, the 
limiting cases of l/TiT [Eq. (fT3|) ] are given by 



(Tir)o 

TiT 



{TiT)o 
T-iT 



13 f^T.) 



TT^ \ a 
for T^Tx 

TT^ \ a 
for T > Tx. 



T/Tx 



(15a) 



(15b) 



The NMR relaxation rate per unit temperature cal- 
culated from the spin fluctuation model is a mono- 
tonic function of temperature. In the high-temperature 
regime l/TiT has a T^^ dependence while in the low- 
temperature regime it is linear in T with a negative slope. 
Thus one realizes immediately that the data for temper- 
atures below Jnmr is not consistent with the predictions 
of the spin fluctuation theory as it has a positive slope. 
We will return to discuss this regime latter. We begin by 
investigating the high temperature regime, T > Tnmr- 

We fit the 1/TiT expression for T :s> lfTtb|) 
to the experimental data of De Sotoi^ for k- 
(ET)2Cu[N(CN)2]Br between Tnmr and 300 K with 
(l/Tir)o, /3[C(r^)/a]^, and as free parameters. It is 



not possible to obtain (3 and (,(Tx)/a unambiguously from 
fitting to 1/TiT data because the model depends sensi- 
tively only on the product (}[^{Tx) / a]'^ (see Eq. (|15b|) ). 
The parameters from the fits are tabulated in Table[T]and 
the results are plotted in Fig. [TJ The use of Eq. (jl5bp to 
fit l/TiT data is justified post hoc since is found to be 
2 — 6 times smaller than Tnmr- We have also checked this 
by plotting the full theory (without taking the T ^ 
limit) for K-(ET)2Cu[N(CN)2]Br in Fig. [Th, where there 
is Korringa ratio data (see Fig. [J) and thus we can de- 
termine P and ^{Tx)/a individually. It can be seen from 
Fig. Wp that the disagreement between the full theory 
and the high temperature approximation is smaller than 
the thickness of the curves, therefore this approximation 
is well justified. It will be shown in Section II B that the 
correlation length is indeed large thus providing further 
justification for the use of Eq. (jl5bp here. 

The model produces a reasonably good fit to the 
experimental data on K-(ET)2Cu[N(CN)2]Br^^ between 
Tnmr, the temperature at which 1/TiT is maximum, 
and room temperature. In the high temperature regime 
(e.g., around room temperature), l/TiT has a very weak 
temperature dependence, indicating weakly correlated 
spins. The large enhancement of 1/TiT can be under- 
stood in terms of the growth of the spin fluctuations: as 
the system cools down, the spin-spin correlations grow 
stronger which allows the nuclear spins to relax faster 
by transferring energy to the rest of the spin degrees of 
freedom via these spin fluctuations. Strong spin fluc- 
tuations, measured by large values of P[^{Tx) / a]"^ , are 
not only present in K-(ET)2Cu[N(CN)2]Br but also ob- 
served in other materials such as the fully deuterated 
K-(ET)2Cu[N(CN)2]Br {which wiU be denoted by K(d8)- 
(ET)2Cu[N(CN)2]Br} and K-(ET)2Cu(NCS)2. The re- 
sults of the fits for K(d8)-(ET)2Cu[N(CN)2]Br and k- 
(ET)2Cu(NCS)2 are shown in Fig. [TJ The parameters 
that produce the best fits are also tabulated in Table [B 
This suggests that strong spin fluctuations are present 
in these charge transfer salts. In all cases studied here, 
strong spin fluctuations are evident from the large value 
of/3K(T,)/a]2. 

The nature of the spin fluctuations, i.e., whether they 
are antiferromagnetic or ferromagnetic, cannot, even in 
principle, be determined from the analysis on 1/TiT. 
Both cases yield the same l/TiT [see Eq. ^ and Sec 
II. C. 2] because the nuclear spin relaxation rate is ob- 
tained by summing all wave vector contribution in the 
first Brillouin zone. However, in the next section we will 
use the Korringa ratio to show that the spin fluctuations 
are antiferromagnetic. 

Below Tnmr, the calculated 1/TiT continues to rise 
while the experimental data shows a decrease in the nu- 
clear spin relaxation rate per unit temperature. However, 
the data does not reach a constant 1/TiT as expected 
for a Fermi liquid. This indicates that the physics below 
Tnmr is dominated by some other mechanism not cap- 
tured by the spin fluctuation, Fermi liquid, or DMFT 
theories. One possibility is a pseudogap opens up at 
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Material 


Ref. 


(l/TiT)o 


(K) 




K-Br 


MayafTre [18] 


0.09 ± 0.01 


6.5 ±5.5 


290 ±250 


K-Br 


De Soto [19] 


0.02 ± 0.01 


20 ±10 


680 ±430 


d8-Br 


Miyagawa [20] 


0.04 ±0.01 


6.2 ±3.5 


85 ±65 


K-NCS 


Kawamoto [21] 


0.06 ±0.01 


11 ±2.6 


110 ±89 



TABLE I: The parameters obtained from the fits which are 
used to produce Fig. [1] Evidence for strong spin fluctuations 
come from the large value of l3[S,{Tx) / which are present for 
all the materials tabulated above. In the table K-Br, d8-Br, 
and K-NCS are used as shorthand for K-(ET)2Cu[N(CN)2]Br, 
K(d8)-(ET)2Cu[N(CN)2]Br, and k-(ET)2Cu(NCS)2 respec- 
tively. 



Tnmr which suppresses the DOS at the Fermi level. Since 
l/TiT - ~P^{Ef) [c.f., Eq. ((Xl6|) ]. a decrease in DOS will 
naturally lead to the suppression of \/TiT. One might 
argue that the discrepancy between the theory and ex- 
periments below Tnmr stems from our assumption of a 
q-independent hyperfine coupling in the \/TiT expres- 
sion. However, in section IIII Dl we will show that the 
Knight shift is also inconsistent with the predictions of 
the spin fluctuation model below Tnmr- While including 
the appropriate q-dependence of the hyperfine coupling 
might change the temperature dependence of 1/TiT, it 
certainly cannot affect the temperature dependence of 
the Knight shift as can be seen from Eq. 



B. The Korringa Ratio 

In the previous section we compared the prediction of 
the spin fluctuation model for 1/TiT to the experimen- 
tal data and obtained good agreement with the data be- 
tween Tnmr and 300 K. However, we were not able to 
determine /3 and ^{Tx)/a unambiguously because 1/TiT 
is sensitive only to the product (3[^{Tx) / a]'^ . We were 
also not able to determine whether antiferromagnetic or 
ferromagnetic spin fluctuations are dominant. We re- 
solve these by studying the Korringa ratio /C. It has 
previously been pointed out that antiferromagnetic (fer- 
romagnetic) fluctuations produce a Korringa ratio that 
is larger (less) than onei^i^ We have also shown in Sec- 
tion II that in the limit of large correlation length, the 
Korringa ratio behaves like (^/a)^ > 1 for antiferromag- 
netic spin fluctuations and like (a/^)^ < 1 for ferro- 
magnetic spin fluctuations. The Korringa ratio data for 
K-(ET)2Cu[N(CN)2]Br (see Fig [2]) is significantly larger 
than one at all temperatures which shows that antiferro- 
magnetic fluctuations dominate. With this in mind, we 
study the antiferromagnetic spin fluctuation model. 

First we note that Kg, given by Eq. has a 

weak temperature dependence because C = 2{£^{Tx) / a)"^ 
is generally larger than unity and T/T^. In the limit 
of large correlation lengths, the second term inside the 
square bracket in the expression for Ks given in Eq. ([Te 



can be approximated as {1 + C + T/Tx) "'^ ~ C ^ and the 
Knight shift wiU be given by ~ (i^a)o(l + VP/{27r^)) 
which is temperature independent. We use this temper- 
ature independent Knight shift to calculate the Korringa 
ratio /C 



1 



(16) 



/C, 



where the prefactor ICq is given by Eq. p4p . 

We flt Eq. (fTB]) to the experimental Korringa data for 
K-(ET)2Cu[N(CN)2]Br.i9 The result is plotted in Fig. H 
In this expression we have three parameters, /3[^(Ta;)/a]^, 
Tx, and ^/P, two of which, /3[^(Ta;)/a]^ and Tx, have been 
determined from fitting l/TiT. There is only one remain- 
ing free parameter in the model, which can then be 
determined unambiguously from the Korringa fit which 
yields /3 = 60 ± 20. This value of /? implies that the an- 
tiferromagnetic correlation length £,{T) — 3.5 ± 2.5a (a 
is the unit of one lattice constant) at T = 50 K. This 
value is in the same order of magnitude as the value of 
the correlation length estimated in the cuprates.''^ 

The Korringa ratio data are well reproduced by the 
antiferromagnetic spin fiuctuation model when T > 
Tnmr- This is again consistent with our earlier con- 
clusion that the spin fluctuations have antiferromagnetic 
correlations. A large Korringa ratio^ii^^ has previously 
been observed in the cuprates indicating similar anti- 
ferromagnetic fluctuations in these systems. The Ko- 
rringa ratio has also been measured in a number of 
heavy fermion compounda^^i^liZ^ Similar antiferromag- 
netic fluctuations, like those observed in the cuprates 
and organics, are present in CeCu2Si2 and. The Kor- 
ringa ratio of this material has a value of 4.6 at 100 mK 
(Ref. [73). In contrast, YbRhzSiz^ and CeRuaSi^I^, 
show strong ferromagnetic spin fluctuations as is evident 
from the Korringa ratio less than unity. In Sr2Ru04^ 
the Korringa ratio is approximately 1.5 at 1.4 K. Upon 
doping with Ca to form Sr2_a:Caa;2R.u04, the Korringa 
ratio becomes less than one which indicates that there 
is a subtle competition between antiferromagnetic and 
ferromagnetic fluctuations in these ruthenates. 



C. The Antiferromagnetic Correlation Length 

It is important to realize that the spin fluctuation for- 
malism can be used to extract quantitative information 
about the spin correlations from NMR data. For exam- 
ple, the fits presented in Figs. [1] and [5] allow us to esti- 
mate the antiferromagnetic correlation length. From the 
fit for K-(ET)2Cu[N(CN)2]Br (Table[I| we found that the 
antiferromagnetic correlation length £^{T)/a = 3.5 ± 2.5 
at r = Tnmr = 50 K. In order to understand the phys- 
ical significance of this value of ^(T) it is informative 
to compare this value with the correlation length for 
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the square-^ and triangular— lattice antiferromagnetic 
Hciscnberg models. 

It has been shown^*' that, on the square lattice, the 
antiferromagnetic Heisenberg model has a correlation 
length of order ^{T)/a ~ 1 for T = J and of order 
^(r)/a - 30 for T = 0.3 J. On the other hand for the 
antiferromagnetic Heisenberg model on the isotropic tri- 
angular lattice, the correlation length is only of order 
a lattice constant at T = Thus the correlation 

length, ^{T)/a = 3.5 ± 2.5, obtained from the analysis 
of the data for K-(ET)2Cu[N(CN)2]Br is reasonable and 
places the materials between the square and isotropic an- 
tiferromagnetic Heisenberg model as has been argued on 
the basis of electronic structure calculations4i^>22 

One of the best ways to measure antiferromagnetic cor- 
relation length is by inelastic neutron scattering experi- 
ments. To perform this experiment, one needs high qual- 
ity single crystal. Unfortunately, it is difficult to grow 
sufficiently large single crystals for k-(ET)2^; however, 
recently some significant progress has been madei^ An- 
other way to probe the correlation length is through the 
spin echo experiment. The spin echo decay rate I/T2 is 
proportional to the temperature dependence correlation 
length [see Eq. ([9])] so measurements of T2 would give us 
direct knowledge on the nature of the correlation length. 
To the authors' knowledge there is no spin echo decay 
rate measurement on the metallic phase of the layered 
organic materials at the present time thus it is very de- 
sirable to have experimental data on T2 measurement to 
compare with the value of £,{T) we have extracted above. 



D. The Knight Shift 

As we pointed out in Section 11 the Knight shift Kg will 
generally have a weak temperature dependence through- 
out the whole temperature range and so, thus far, we 
have neglected its temperature dependence. However, it 
is apparent from Eq. (jl3[) that for any choice of parame- 
ter values {(3,£,{Tx)/a, and T^}, Kg will always increase 
monotically as the temperature decreases. Therefore the 
temperature dependence of the Knight shift potentially 
provides an important check on the validity of the spin 
fluctuation model. However, in the following discussion 
one should recall the caveats (discussed in section IH Ap 
on the validity of the calculation of the Knight shift stem- 
ming from the assumption that the dynamics of the long 
wavelength part of dynamical susceptibility relax in the 
same manner as a Fermi liquid. 

In contrast to the prediction of the spin fluctuation 
model the experimental data, e.g., those measuredi^ on 
K-(ET)2Cu[N(CN)2]Br (reproduced in Fig. g]), show that 
Kg decreases slowly with decreasing temperature which 
then undergoes a large suppression around Tk^ ~ 50 K. 
It should be emphasized here that Tk^ is approximately 
the same as Tnmr, the temperature at which 1/TiT is 
maximum. 

Since it is not possible to explain any of the NMR 



data below Tnmr in terms of the spin fluctuation model 
within the approximations discussed thus far we focus on 
the temperature range between 50 K to 300 K just as we 
did for the analysis of 1/TiT. Even in this temperature 
range, there is a puzzling discrepancy between theory and 
experiment: the experimental data decreases slowly with 
decreasing temperature while the theoretical calculation 
predicts the opposite. We will argue below that this dis- 
crepancy arises because the data are obtained at con- 
stant pressure while the theoretical prediction assumes 
constant volume. Since the organic charge transfer salts 
are particularly soft, thermal expansion of the unit cell 
may produce a sizeable effect to the Knight shift and may 
not be neglected. 

Following Wzietek et al.^ we attempt to make an 
estimate on the correction of the Knight shift for n- 
(ET)2Cu[N(CN)2]Br due to thermal expansion. Let us 
define {T, V) as the constant volume spin susceptibil- 
ity as a function of temperature. The measured sus- 
ceptibility is then given by a constant pressure suscep- 
tibility = X^s[T,V{T,P)] while the theoretical sus- 
ceptibility is given by a constant volume susceptibility 
= x'"s[T,V{T = 0,P)]. The correction to the spin 
susceptibility is then given by 



Ax 



dT' 



Jo 



fVdP\ / dV 
'dPjj,, \ dV J J,, \VdT' 



(17) 



The Knight shift is directly proportional to the spin sus- 
ceptibility, Eq. (Pc)) . which allows us to write the correc- 
tion to the Knight shift as 



•T 



(18) 



dV 
VdT' 



where is the (experimentally obtained) isobaric 
Knight shift, is the (calculated) constant volume 
Knight shift, {VdP/dV)T is the isothermal compress- 
ibility, and {dV/VdT)p is the linear thermal expansion. 
It is hard to obtain an accurate estimate for AJCj be- 
cause there are no complete sets of data for iff, isother- 
mal compressibility, and thermal expansion as a function 
of temperature and pressure for the «;-(ET)2-X family. 
However a rough estimate for Aifs may be made using 
the available experimental data. 
In Appendix [B] we estimate that 



-3X10- bar- 



V 



and 



dP_ 
dV 

Vdf 



-10^ bar, 



4 iz-l 



IQ-'' K 



Combining these order of magnitude estimates we are 
able to obtain a rough estimate on AKg which can be 
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FIG. 4; The temperature dependence of the (constant pressure) Knight shift as measured by K-(ET)2Cu[N(CN)2]Br by De 
Soto et al.^^ (filled squares) and the corrected Knight shift obtained by taking into account thermal expansion of the lattice 
(half-filled diamonds), i.e., the constant volume Knight shift. The temperature at which the Knight shift decreases rapidly is 
about the same temperature at which 1/TiT is suppressed (see Fig. [T|, i.e. Tk^ ~ Tnmr. This suggests that 1/T\T and Ks are 
suppressed by the same physics. In the limit of large correlation lengths, the spin fluctuation model, Eq. ((13} , predicts a slowly 
varying Knight shift which is almost temperature independent (solid line). The discrepancy between theory and experiment 
arises because the model calculates constant volume Knight shift while the experiment measures constant pressure Knight 
shiftjSi Since k-(ET)2X is soft there will be a sizeable effect on Kg due to the large thermal expansion. The correction to the 
experimental Knight shift, by taking into account these effects, was calculated by using Eq. (|18|l (half filled squares). However, 
we stress that the lack of compressive measurements of pressure and temperature dependence of the Knight shift, isothermal 
compressibility, and thermal expansion means that this correction is no better than an order of magnitude estimate. However, 
our estimate indicates that there correction is large enough that the data above Tk^ cannot be shown to be in disagreement 
with the spin fluctuation model. Below Tk^ there is a clear disagreement between the theory and data, in section |V] we argue 
that this is because a pseudogap opens at Tnmr. In the key to this figure «:-(ET)2Cu[N(CN)2]Br is abbreviated as K-Br. 




written as ~ - 0.3T for T in Kelvin. The re- 
sult is plotted in Fig. 21 It is clear from the figure that 
our rough estimate has already produced a non trivial 
correction to the Knight shift. The Knight shift changes 
from having a positive slope in the raw data to exhibit 
a rather small negative slope between Tk^ ~ 50 K and 
room temperature when the corrections to account for 
the thermal expansion are included. The correction be- 
comes small below about 50 K. The lattice expansion 
clearly has a significant effect on the measured Knight 
shift. To remove this effect one would need to either 
measure the Knight shift at constant volume or pursue 
an experiment in which the pressure dependence of Ks^ 
isothermal compressibility, and thermal expansion [c.f. 
Eq. are measured simultaneously to accurately de- 

termine AfCj. Given the large uncertainty in AfCj we 
take Kg to be constant for temperatures above 50 K in 
the rest of this paper. This is clearly the simplest as- 
sumption, it is not (yet) contradicted by experimental 
data, and, perhaps most important, any temperature de- 
pendence in the Knight shift is significantly smaller than 
the temperature dependence of 1/TiT. 

Regardless of the valuse of AKg, the Knight shift cal- 
culated from the spin fluctuation model is inconsistent 
with the experimental data below Tr-^ ~ 50 K (see Fig. 
2]) . The calculated Kg shows a weakly increasing Kg with 
decreasing temperature, while the measured Kg is heavily 



suppressed below 50 K. One important point to empha- 
size here is that the temperature dependence of Kg will 
not change even if one uses the fully q-dependent A{q) 
since Kg [see Eq. (Pel) ] only probes the q = compo- 
nent of the hyperfine coupling and susceptibility. Thus, 
putting an appropriate q-dependent hyperfine coupling 
will not change the result for Kg (although it might give 
a better description for l/TiT). This provides a com- 
pelling clue that some non-trivial mechanism is respon- 
sible to the suppression of l/TiT, Kg, and /C below 50 
K. 

We have not addressed how the nuclear spin relaxation 
rate is modified by the thermal expansion of the lattice. 
Since the organic compound is soft, it is interesting to 
ask if there is a sizeable effect to 1/TiT. Wzietek et 
al§^ have performed this analysis on quasi- ID organic 
compounds whose relaxation rate in found to scale like 
Xs- One can straightforwardly derive the effect of volume 
changes from the Hubbard model. If one uses the relation 
l/TiT - xl and assumes fixed U and t, then l/TiT - 
l/V"^ will follow. However, it is clear from the phase 
diagram of the organic charge transfer salts (Fig. O and 
Ref. Q) that there is a rather large change in U and t 
for even small pressure variations. Therefore, there is 
no obvious relationship between 1/TiT and Xs for the 
quasi-2D organics and it is not clear how the imaginary 
part of the susceptibility x"{^.^^)^ which enters l/TiT, 
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is effected by thermal expansion and lattice isothermal 
compressibility. More detailed experiments are clearly 
needed to determine the effect of thermal expansion of 
the lattice on the measured relaxation rate. 



IV. SPIN FLUCTUATIONS IN THE MOTT 
INSULATING PHASE OF k-(ET)2CU2(CN)3 

Recent experiments on k-(ET)2Cu2(CN)3 by Shimizu 
and collaborators2M^i^i22 have generated a lot of 
interest This is because the Mott insulat- 

ing phase of this material appears to have a spin liquid 
ground state, that is a state which does not have mag- 
netic ordering (or break any other symmetry of the nor- 
mal state) even though well- formed local moments exist. 
This is very different from the Mott insulating phases of 
the other n salts, such as k-(ET)2Cu[N(CN)2]C1, which 
clearly shows antiferromagnetic ordering^^ at low tem- 
perature and ambient pressure. An elegant demonstra- 
tion of these two different ground states is provided by 
susceptibility measurements!^ the susceptibility of k- 
(ET)2Cu[N(CN)2]C1 exhibits an abrupt increase around 
25 K which marks the onset of Neel ordering while the 
susceptibility of k-(ET)2Cu2(CN)3 shows no sign of a 
magnetic transition. The transition to a magnetically or- 
dered ground state reahzed in k-(ET)2Cu[N(CN)2]C1 is 
also demonstrated by the splitting of NMR spectra 
below the transition temperature^^ The difference in 
the ground states of k-(ET)2Cu[N(CN)2]C1 and k- 
(ET)2Cu2(CN)3 appears to be connected with the fact 
that there is significantly greater frustration in k- 
(ET)2Cu2(CN)3 (for which t' /t - 1) than there is in k- 
(ET)2Cu[N(CN)2]Cl (for which i'/t- 0.7). Geometrical 
frustration alone is not sufficient to explain the absence 
of magnetic order in k-(ET)2Cu2(CN)3 because a Heisen- 
berg model on an isotropic triangular lattice is known 
to exhibit a magnetically ordered ground state, i.e. 120° 
state. It may be that the proximity to the Mott transition 
plays an important role in allowing the absence of mag- 
netic ordering at low temperatures in k-(ET)2Cu2(CN)3. 
One possible explanation for the existence of a spin liq- 
uid ground state is there are ring exchange terms in the 
Hamiltonian arising from charge fluctuations which has 
been studied by several groups i^i^SiSi 

The NMR relaxation rate in k-(ET)2Cu2(CN)3 (Ref. 
Issi and Fig. [5]) shows a similar temperature dependence 
to that in, for example, K-(ET)2Cu[N(CN)2]Br. l/TxT is 
enhanced over the Korringa-like behavior with a peak at 
Tnmr ~ 10 K below which it exhibits a large decrease. 
However the Knight shift Kg in k-(ET)2Cu2(CN)3 is 
quite different to that in K-(ET)2Cu[N(CN)2]Br (com- 
pare Figs. HandE]). In k-(ET)2Cu2(CN)3, Kg increases 
as the temperature is lowered from room temperature 
until it reaches a broad maximum around Tk^ ~ 30 — 50 
K below which it drops rapidly. In contrast. Kg in k- 
(ET)2Cu[N(CN)2]Br shows a weak temperature depen- 
dence down to Tk, below which it undergoes a sharp de- 



crease (see Fig. |4]). Another difference is Tnmr is consid- 
erably lower than Tk^ in k-(ET)2Cu2(CN)3 whereas they 
are roughly the same in K-(ET)2Cu[N(CN)2]Br. This 
suggests that whereas the suppression of 1/TiT below 
Tnmr and Ks below Tk, in K-(ET)2Cu[N(CN)2]Br prob- 
ably has a common origin; in k-(ET)2Cu2(CN)3 the ori- 
gin of the suppression of l/TiT below Tnmr is different 
from the origin of the broad maximum in Ks at Tr-^. 
Note that the fact that JC > 1 shows that the spin fluc- 
tuations are antiferromagnetic, this is rather interesting 
given the importance of Nagaoka ferromagnetism on the 
triangular latticei^ 

Given the reasonable agreement between the antiferro- 
magnetic spin fluctuation model with the NMR data on 
K-(ET)2Cu[N(CN)2]Br (down to Tnmr 50 K), we ap- 
ply the same formalism to k-(ET)2Cu2(CN)3. A slight 
modification to the spin fluctuation model is necessary 
since, unlike the other k salts studied in this paper, k- 
(ET)2Cu2(CN)3 is an insulator. Therefore we clearly 
cannot use the Fermi liquid form of xlw- The sim- 
plest approximation is that the dynamic susceptibility 
given in Eq. ([3]) will only consist of XAF(q, 'j-')- In 
the region where £_{T)/a is large, XQ = and 
wgp = Q;'(^/a)~^ where a and a' are temperature inde- 
pendent constants and a is the lattice spacing. Within 
these approximations the nuclear spin relaxation rate and 
Knight shift are given by 



1 



TiT/o l + (Qa)2(e/a)2 

(e/a)'-" 



with 



{Ks)a 
iKs)o 



l + {Qa)^^/ay 



2afcs|A|2 



(19) 



(20) 



where Q = {Q,Q) is the finite wave vector on which we 
assume the susceptibility to peak. Again, we take the 
temperature dependence of the correlation length to be 

Following the same approximation scheme as before 
(outlined in Section II B), we assume a relaxational dy- 
namics of the spin fluctuations, which are described by 
a dynamic critical exponent z = 2, and a mean field 
critical exponent 77 = 0. Within these approximations 
wgp = a(a/^)^ and XQ = Q;'(C/'^)^- The nuclear spin 
relaxation rate and Knight shift are then given by 



1 

Ks 



(l/TiT)oC2 



(T/T, + 1)2 + {QayCiT/T, + 1) 
{Ks)qC 



1 -f {QafC -t- T/T, 



(21) 



We work in the high temperature approximation for 
1/TiT - following a procedure similar to that em- 
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FIG. 5: [Color online] Comparison of the spin fluctuation theory with the measureti^ temperature dependence of the nuclear 
spin relaxation rate per unit temperature, l/TiT (left panel), Knight shift, Ks (center panel), and Korringa ratio (right panel) 
of the Mott insulating phase of k-(ET)2Cu2(CN)3 [abbreviated as k-(CN)3 in the keys to the figures]. The spin fluctuation 
model is in good agrement with the measured 1/TiT, but does not describe the Knight shift (and hence IC) well. We have 
also checked that using the Fermi liquid for Xlw does not improve the fit to the Ks data, and this fit is shown as a dashed 
line. This fit to the data is clearly worse than simply setting xlw ~ 0. Also note that the peak in l/TiT, Tnmr, is at a lower 
temperature than the maximum in the Knight shift, Tk,- This behavior is qualitatively different from that of the other k salts 
where Tnmr ~ Tk, (see Figs. [T]and|3|. This suggests that in «-(ET)2Cu2(CN)3 the origin of the 1/T\T suppression is different 
from physics that gives rise to the maximum in Ks . The parameter values from the lines of best fit shown in this figure are 
reported in Table|ITl The fact that the model gives a reasonably good fit to 1/TiT but not to the Ks data suggests that the spin 
fluctuation model fails to correctly account for the long wavelength physics, but suggests that the model correctly describes the 
physics around a peak in x(q,<^) which dominates the integral over the first Brillouin zone and thus 1/TiT [c.f., Eq. (I2a|) ]. This 
is rather surprising as the correlation length is less than one lattice constant at T = 50 K. This result is clearly inconsistent 
with the initial assumption that the long wavelength susceptibility is dominated by a peak in the dynamic susceptibility at 
a finite wave vector. Finally we note that the fact that /C > 1 shows that the spin fluctuations are antiferromagnetic, this is 
rather interesting given the importance of Nagaoka ferromagnetism on the triangular lattice.— 



ployed to obtain Eq. (jl5bp . If the dynamic sus- 
ceptibility is strongly peaked at q — Q then the pa- 
rameter 2{Qa)'^{£^{Tx)/a)'^ is much larger than 1 so 
{QafC{T/T^ + 1) is larger than [T/T.;^ + 1)^ and we 
can just keep the term proportional to (Qa)^ in the de- 
nominator which allows us to write 1/TiT as 

J_ (l/Tir)o[V2C(T.)/a]^ 

TiT" (Qa)2(r/r, + 1) ^ ^ 

We use the expressions for l/TiT given in Eq. ([22|) and 
for Ks in Eq. ^ to fit the k-(ET)2Cu2(CN)3 data.55 
We assume the susceptibility has a strong peak at Q 
27r/3.^^^ The parameters of the best fit are reported in 
Table [II] and the results are plotted in Fig O While the 
spin fluctuation model produces a reasonably good fit to 
the 1/TxT above Tnmr ~ 10 K, it does not reproduce 
Ks data well as can be seen from the upward curvature 
in the fit in contrast to the data which shows a slight 
downward curvature. We also performed the fit to k- 
(ET)2Cu2(CN)3 data using the most general forms Eq. 

and taking z and 77 as free parameters. Good fits to 
l/TiT and Kg can be obtained with 77 '--^ 1 and z ^ 2 but 
this gives us so many free parameters that the value of 
such fits must be questioned. 

An important fact is that Ks probes the long wave- 
length dynamics. We have set xvn{'^) fo zero in order 
to make the simplest possible assumption about the in- 
sulating state. The data indicate that this assumption 
is probably incorrect. Recently Zheng et al^ used a 
high temperature series expansion to calculate the uni- 



Parameter 


Fit results 


(l/TiT),, (s-^K-i) 


220 ±11 


{Ks)o 


8100 ± 720 


(K) 


40 ± 4 




0.3 ± 0.1 



TABLE II: The parameters obtained from the best fits to 
l/TiT and Ks data in k-(ET)2Cu2(CN)3. These parameters 
are used to produce Fig. [5] The antiferromagnetic correlation 
length £,{Tx)/a is short ranged consistent with the significant 
frustratio n^'^ '^^'^^ present in this material. 



form spin susceptibility (which is the same as the Knight 
shift apart from a constant of proportionality) for the 
Heisenberg model on a triangular lattice, applied it to k- 
(ET)2Cu2(CN)3. They obtained a good agreement with 
the experimental data. The spin fluctuation model de- 
scribed here can be viewed as a different route to under- 
stand the same experiment. The discrepancy between the 
spin fluctuation model and the data suggests a failure of 
our implicit assumption that the long wavelength physics, 
which determines Ks , is dominated by a peak in the dy- 
namic susceptibility due to spin fluctuations. This is con- 
sistent with the fact that we find that ^(T) ~ 0.2—0.4 lat- 
tice spacings at T = 50 K which clearly disagrees with our 
initial assumption that ^(T) ^ a. This begs the ques- 
tion what physics dominates the long wavelength physics 
both in the series expansions and in the real material? 
The low temperature properties of k- 
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(ET)2Cu2(CN)3 are clearly inconsistent with a magnetic 
ordered ground state. For a two-dimensional quantum 
spin system with an ordered ground state, the low tem- 
perature properties are captured by the non-linear sigma 
model. The observed temperature dependence of l/Ti 
and the spin echo rate I/T2 foUow 1/Ti cx r^/2^(T), 
and I/T2 cx T^^{T), where the correlation length ^(T) is 
given by'''^ 



= 0.021 



T 



1/2 



cxp 



T 



(23) 



where c is the spin wave velocity and ps is the spin stiff- 
ness. In the quantum critical regime/^ l/TiT ^ yj^i 
and I/T2 ~ T''-! [c.f., Eq. ^] where 77 is the anoma- 
lous critical exponent associated with the spin-spin cor- 
relation function whose value is generally less than 1. 
Thus for a magnetically ordered state, which can be 
well described by 0(N) non linear sigma model, both 
1/TiT and I/T2 should increase with decreasing tem- 
perature. For K-(ET)2Cu2(CN)3 Shimizu et al.^^ found 
that l/TiT ~ and l/r2 - constant from 1 K down 
to 20 mK which suggests the critical exponent rj > 1. 
The nuclear spin relaxation rate decreases with decreas- 
ing temperature. Such a large value of z is what occurs 
for deconfined spinous 



V. UNCONVENTIONAL COHERENT 
TRANSPORT REGIME AND THE 
K-(ET)2X PHASE DIAGRAM 



In Section Hi Bl we discussed the DMFT description of 
the crossover from a bad metal to a Fermi liquid. DMFT 
successfully predicts the unconventional behaviors ob- 
served in a number of experiments on the k-{FiT)2X salts. 
These include the resistivity, thermopower, and ultra- 
sound velocity. The unconventional behaviors seen in 
these measurements are associated with the crossover 
from bad metallic regime to a renormalized Fermi liquid 
in the DMFT picture. While DMFT gives reUable pre- 
dictions for the transport properties r^i^^i^i^^ it is not 
able to explain the loss of DOS observed in the nuclear 
spin relaxation rate and Knight shift. Thus the NMR 
data suggest that the coherent transport regime is not 
simply a Fermi liquid, contrary to what has previously 
been thought. 

To illustrate the nature of the low temperature param- 
agnetic metallic state more qualitatively, it is instructive 
to study how the bad metal-coherent transport crossover 
is related to the loss of DOS. Therefore we have investi- 
gated the relationship between Tp^x2, the temperature 
at which the resistivity deviates from behavior; T^„/^, 
the temperature at which a dip in the ultrasonic velocity 
is observed; and Tnmr, the temperature at which 1/TiT 
(and Ks) is maximum which appears to mark the onset 
of a loss of DOS. In Figure [5] we plot Tp^r2, T^^/^,, and 
Tnmr as measured by several different groups for various 
salts against Tc which serves well as a single parameter to 



Material 


Tc{P) 




Tav/v{P) 


rNMR(-P) 


K-C\ 


[38^ 


[25] 


[38] 


[95] 


K-Br 


[37^ 


M 


[37] 


[18] 


K-NCS 


[98j 




[37] 




K-(CN)3 


[86j 


[86] 




[85] 



TABLE III: The references from which the pressure depen- 
dence of different temperature scales for different materials 
used to produce Fig. [S]were taken. The notation is the same 
as that given in Fig. [6] Tp^j-i is the temperature at which 
the resistivity deviates from a behavior, T^v/v is the tem- 
perature at which a dip in the ultrasound velocity is observed, 
and Tnmr is the temperature on which 1/T\T peaks. In the 
table K-Cl, K-Br, and k-NCS and fi:-(CN)3 are used as short- 
hand for fi:-(ET)2Cu[N(CN)2]Cl, fi:-(ET)2Cu[N(CN)2]Br, k- 
(ET)2Cu(NCS)2, and k-(ET)2Cu2(CN)3 respectively. 



characterize both the hydrostatic pressure and the varia- 
tion in chemistry (or 'chemical pressure'). This analysis 
is complicated by the necessity of comparing pressures 
from different experiments. Our procedure for dealing 
with this issue is outlined in the Appendix [Cl and we 
stress that the large error bars in Fig. [5] are due to the 
difficulties in accurately measuring pressure rather than 
problems in determining Tc, Tp^^T^, Ja^/^, or Tnmr- 

It is clear from Fig. [6] that the data for k- 
(ET)2Cu[N(CN)2]Cl, K-(ET)2Cu[N(CN)2]Br, and k- 
(ET)2Cu(NCS)2 fall roughly onto a single curve. This 
suggests that Tnmr coincides with Tp^x^ and T^y/^. 
Thus the loss of DOS, associated with Tnmr, occurs 
around the temperature at which the crossover from bad 
metal to coherent transport regime takes place. The 
loss of DOS observed in l/TiT and Kg is not what one 
would expect for a Fermi liquid; therefore the coherent 
intralayer transport regime is more complicated than a 
renormalized Fermi liquid. This must result from non- 
local correlations which are not captured by DMFT since 
DMFT captures local correlations exactly. One possible 
explanation for the loss of DOS is the opening of a pseu- 
dogap. 

Another important point to emphasize from Fig. 
[6] is the appearance of a second trend formed by 
K-(ET)2Cu2(CN)3 which is clearly distinct from the 
trend of the data points from k-(ET)2Cu[N(CN)2]C1, 
K-(ET)2Cu[N(CN)2]Br, and k-(ET)2Cu(NCS)2. This 
shows that the spin fluctuations in k-(ET)2Cu2(CN)3 are 
qualitatively different from those in the other k- 
{KT)2X salts. Of course, qualitative differences are not 
entirely unexpected due to the spin liquid rather than 
antiferromagnetic ground state in k-(ET)2Cu2(CN)3. It 
has recently been argued that the differences in the 
spin fluctuations in k-(ET)2Cu2(CN)3 will lead this 
material to display a superconductivity with a differ- 
ent symmetry of the order parameter than the other 
K-(ET)2X saltsi^i^ This result shows that the spin 
fluctuations are indeed qualitatively different in k- 
(ET)2Cu2(CN)3. 
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FIG. 6: [Color online] The relationship between different temperature scales for a range of organic charge transfer salts. The 
superconducting transition temperature Tc is used to parameterize the proximity of the material to the Mott transition (Tc 
decreases as one moves further away from the Mott transition). A plot of ,T^^/y, and Tnmr against Tc for several 

K-(ET)2X salts shows that the peak in 1/T\T, Tnmr, occurs at the same temperature as the crossover form a bad metal to 
the coherent transport regime measured in transport {Tp^j,2) and ultrasonic attenuation (Tak/i,) experiments. is the 

temperature at which the resistivity deviates from a behavior, Tad/d is the temperature at which a dip in the ultrasound 
velocity is observed. The left panel shows the data for the k-(ET)2X family in the metallic phase {fi:-(ET)2Cu[N(CN)2]Cl and 
K-(ET)2Cu2(CN)3 under pressure, K-(ET)2Cu[N(CN)2]Br and k-(ET)2Cu(NCS)2} while the right panel shows the data for the 
K-(ET)2X family in the insulating phase [k-(ET)2Cu[N(CN)2]C1 and «;-(ET)2Cu2(CN)3 at ambient pressure]. In the metallic 
phase we use Tc as a single parameter to characterize the effect of chemical substitution and hydrostatic pressure. This works 
surprisingly weU and the data for k-(ET)2Cu[N(CN)2]C1, K-(ET)2Cu[N(CN)2]Br, and k-(ET)2Cu(NCS)2 is seen to collapse 
roughly onto a single trend, which suggests that the spin fluctuations in the metallic phases are rather similar. In contrast, 
the data for k-(ET)2Cu2(CN)3 fall onto a separate curve which suggests that there are important differences between the 
spin fluctuations in this material and those in other k phase salts. This is perhaps not so surprising in light of the fact 
that k-(ET)2Cu2(CN)3 has a spin liquid round state in its Mott insulating phase while the materials are close to at Neel 
ordered Mott insulating phase. This plot suggests that the pseudogap opens at the same temperature as the crossover from 
bad metal to coherent transport regime. Whether this is because of a deep link between the crossover and the pseudogap 
or because the lack of coherence in the bad metal destroys the pseudogap remains to be seen. Collectively these data show 
that the coherent transport regime is not simply a renormalised Fermi liquid as has previously been thought. It should be 
emphasized that the large error bars are the result of our estimates of the systematic errors produced by equating pressures in 
different experiments. The procedure to obtain the error bars presented in this plot is discussed in Appendix [C] The symbols 
represent both the material and the experiment as follows: flUed symbols correspond to the data for k-(ET)2Cu[N(CN)2]C1, 
open symbols denote K-(ET)2Cu[N(CN)2]Br, open symbols with black dots denote k-(ET)2Cu(NCS)2, and half flUed symbols 
denote k-(ET)2Cu2(CN)3. The square symbols represent Tp^y2 vs Tc, circles represent Tai,/„ vs Tc, and triangles represent 
Tnmr vs Tc. The references from which the data were collected are given in Table HITl 



On the basis of the above analysis we sketch the phase 
diagram of k-(ET)2X , shown in Fig. [H The pseudogap 
phase shows an interesting set of behaviors. On one hand 
it exhibits a loss of DOS as is evident from 1/TiT and 
Kg. On the other hand, it exhibits coherent intralayer 
transport as is shown by the resistivity behavior^'^; 
it also has long lived quasiparticles and a well defined 
Fermi surface clearly seen from de Haas- van Alphcn and 
Shubnikov-de Haas oscillation experiments. One 
framework in which it may be possible to understand 



both of these sets of behaviors is if there is a fluctuat- 
ing superconducting gapi^S This idea has been applied to 
the cupratesji2£ii£i it would be interesting to see whether 
such an approach gives a good description of k-(ET)2X. 
Another interesting observation is that the measurements 
which see the loss in the DOS probe the spin degrees of 
freedom whereas the evidence for well defined quasiparti- 
cles comes from probes of the charge degrees of freedom. 
This may be suggestive of a 'spin gap' which could result 
from singlet formation as in the RVB picture. 
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To date there have been few experiments studying the 
pressure dependence of 1/TiT or Kg- Therefore it is not 
possible, at present, to determine with great accuracy 
where the pseudogap vanishes. The available NMR ex- 
periment under pressure -"^'^^ suggest that at sufficiently 
high pressures the pseudogap onset temperature is lower 
than the incoherent-coherent crossover temperature and 
pressure eventually suppresses the pseudogap altogether. 
We represent the current uncertainty over where the 
pseudogap is completely suppressed by pressure by draw- 
ing a shaded area with a question mark in the phase di- 
agram. It is plausible that the pseudogap vanishes very 
close, if not at the same pressure, to the point where the 
superconducting gap vanishes. This would be consistent 
with RVB calculations^^-^^ which suggest that the pseu- 
dogap and superconducting gap are proportional to each 
other and so should vanish at about the same pressure. 
However, we should stress that there really is not yet 
sufficient data to determine exactly where the pseudogap 
vanishes and admit that our choice is, perhaps, a little 
provocative. 

Clearly a series of careful experiments are required 
to elucidate when Tnmr tends to zero. Understanding 
where the pseudogap vanishes is an important consider- 
ation in light of the number of theories based on a hid- 
den pseudogap quantum critical point in the cupratesJ^ 
Furthermore, the superconducting state in organic charge 
transfer salts far from the Mott transition are highly un- 
conventional. These low Tc organic charge transfer salts 
have unexpectedly large penetration depths^i^ and are 
not described by BCS theoryi^ The possibility that a 
quantum critical point is associated with the pressure 
where Tc goes to zero invites comparison with the heavy 
fermion material CeCoIn5_^Sn2;^*''' in which a quantum 
critical point seems to be associated with the critical dop- 
ing to suppress superconductivity. Thus low Tc organic 
charge transfer salts appear increasingly crucial for our 
understanding of the organic charge transfer salts. ^ 

An important question to address theoretically is why 
Tnmr might coincide with Tp^^^ and T^u/t,. Of course, 
it may be that the two phenomena are intimately con- 
nected. However, another possibility suggests itself on 
the basis of DMFT and RVB calculations. DMFT 
correctly captures the local physics and it is this lo- 
cal physics that dominates the cross-over from a 'bad- 
metal' to coherent in-plane transport. On the other 
hand RVB does not capture this cross-over (because the 
Mott transition is only dealt with at the Brinkmann-Rice 
level^"^) but does capture some of the non-local physics 
which DMFT neglects. The pseudogap is predicted by 
RVB theory to increase in temperature when pressure is 
lowered . ^^'^^ This rise in the pseudogap temperature is 
predicted to continue until the pressure is lowered all the 
way to the Mott transition, in contrast with the observed 
behavior (c.f.. Figs [3] and [6]). However, we conjecture that 
the RVB physics is 'cut off' by the loss of coherence at 
Tp^y2 and T^„/« thus preventing Tnmr, from exceeding 



VI. CONCLUSIONS 

We have applied a spin fluctuation model to study the 
temperature dependences of the nuclear spin relaxation 
rate and Knight shift in the paramagnetic metallic phases 
of several quasi two-dimensional organic charge transfer 
salts. The large enhancement of 1/TiT between Tnmr 
{- 50 K in K-(ET)2Cu[N(CN)2]Br} and room tempera- 
ture has been shown to be the result of strong antiferro- 
magnetic spin fluctuations. The antiferromagnetic corre- 
lation length is estimated to be 3.5 ± 2.5 lattice spacings 
in K-(ET)2Cu[N(CN)2]Br at T = 50 K. The temperature 
dependence of 1/TiT for T>Tnmr from the spin fluctua- 
tion model is qualitatively similar with the predictions of 
DMFT. The spin fluctuations in k-(ET)2Cu[N(CN)2]C1, 
K-(ET)2Cu[N(CN)2]Br, K(d8)-(ET)2Cu[N(CN)2]Br, and 
K-(ET)2Cu(NCS)2 are found to be remarkably similar 
both qualitatively and quantitatively. Strong spin fluc- 
tuations seem to be manifested in materials close to 
Mott transition. Recent NMR experimentsiS^ on k- 
(ET)2Ag(CN)2-H20, which is situated further away from 
the Mott transition, suggests that the spin fluctuations 
in this materials are not as strong as those in the other 
K salts studied here. 

We have also applied the spin fluctuation formalism 
to the strongly frustrated system k-(ET)2Cu2(CN)3. In 
this compound the measured 1/TiT, which probes the 
entire Brillouin zone, agrees well with the predictions of 
the spin fluctuation model for T > Tnmr ~ 10 K. In con- 
trast the measured Knight shift, which only depends on 
the long wavelength physics, is not well described by the 
spin fluctuation model. This suggests that at least one 
of the assumptions made in the spin fluctuation model: 
(i) z — 2, rj — 0, or (ii) x(q, w) is strongly peaked at 
wave vector q = Q; is violated in k-(ET)2Cu2(CN)3 or 
that the model neglects some important long wavelength 
physics. In light of the recent evidence for a spin liq- 
uid ground state in k-(ET)2Cu2(CN)3, in contrast to the 
antiferromagnetic or 'd-wave' superconducting grounds 
states in k-(ET)2Cu[N(CN)2]C1, K-(ET)2Cu[N(CN)2]Br, 
K(d8)-(ET)2Cu[N(CN)2]Br, and k-(ET)2Cu(NCS)2, and 
the greater degree of frustration in k-(ET)2Cu2(CN)3 it is 
interesting that there are such important qualitative and 
quantitative differences between the spin fluctuations in 
K-(ET)2Cu2(CN)3 and those in the other k phase salts. 

The peak of 1/TiT and the suppression of Kg are 
strongly dependent on pressure: they are systematically 
reduced and completely vanish at high pressure (> 4 
kbar)jiS, at high pressure a Korringa-like temperature de- 
pendences of 1/TiT and are recovered for all temper- 
atures. It is clear that high pressures will suppress both 
the antiferromagnetic spin fluctuations which are dom- 
inant above 50 K and the mechanism (presumably the 
pseudogap) which causes drops in 1/TiT and Kg below 
50 K at ambient pressure. 

The large suppression of 1/TiT and Ks below Tnmr 
observed in all the k salts studied here cannot be ex- 
plained by the M-MMP spin fluctuation model. The 



17 



most plausible mechanism to account for this feature is 
the appearance of a pseudogap which causes the suppres- 
sion of the density of states at the Fermi energy. This is 
because at low temperature 1/TiT and Ks are propor- 
tional p^{Ep) and p{Ep), respectively [c.f., Eq. (??)]. 
Independent evidence for the suppression of density of 
states at the Fermi level comes from the linear coefficient 
of specific heat The electronic specific heat probes 
the density of excitations within fc^T of the Fermi en- 
ergy. Any gap will suppress the density of states near 
the Fermi surface which results in the depression of the 
specific heat coefficient 7. Kanoda^^ compared 7 for sev- 
eral of the K-(ET)2X salts and found that in the region 
close to the Mott transition, 7 is indeed reduced. One 
possible interpretation of this behavior is a pseudogap 
which becomes bigger as one approaches the Mott tran- 
sition. However, other interpretations are also possible, 
in particular one needs to take care to account for the 
coexistence of metallic and insulating phases; this is ex- 
pected as the Mott transition is first order in the organic 
charge transfer salts j^i^Si The existence of a pseudogap 
has also been suggested A-(BEDT-TSF)2GaCl4i"*^ from 
microwave conductivity. The reduction of the real part 
of the conductivity ui from the Drude conductivity Udc 
and the steep upturn in the imaginary part of the conduc- 
tivity (72 may be interpreted in term of preformed pairs 
leading to a pseudogap in this material. 

The experimental evidence from measurements of 
1/TiT, ifjj, and heat capacity all seem to point 
to the existence of a pseudogap below Tnmr in 
(ET)2Cu[N(CN)2]Br and k-(ET)2Cu(NCS)2. Thus a 
phenomenological description which takes into account 
both the spin fluctuations which are important above 
Tnmr and a pseudogap which dominates the physics be- 
low Tnmr would seem to be a reasonable starting point 
to explain the NMR data for the entire temperature 
range (clearly superconductivity must also be included 
for T < Tc). We will pursue this approach in our fu- 
ture work. In particular one would like to answer the 
following questions: how big is the pseudogap and what 
symmetry does it have? Is there any relation between 
the pseudogap and the superconducting gap? The an- 
swer to these questions may help put constraints on the 
microscopic theories. 

Future experiments. There are a number of key ex- 
periments to study the pseudogap. The pressure and 
magnetic field dependences of the nuclear spin relaxation 
rate and Knight shift will be valuable in determining the 
pseudogap phase boundary, estimating the order of mag- 
nitude of the pseudogap, and addressing the issue how 
the pseudogap is related to superconductivity. In the 
cuprates, there have been several investigations of the 
magnetic field dependence of the pseudogap seen in NMR 
experiments. For Bi2Sri.6Lao.4Cu06 the nuclear spin re- 
laxation rate does not change will field up to 43 ^22^ 
However, since T* ^ 200 K, one may require a larger 
field to reduce the pseudogap. Similar results were found 
in YBa2Cu408.^^" However, in YBa2Cu307_5 [see espe- 



cially Fig. 6ofRef. |m| a field of order 10 T is enough to 
start to close the pseudogap. Mitrovic et ali^ interpret 
this observation in terms of the suppression of 'd-wave' 
superconducting fluctuations. 

The interlayer magnetoresistance of the cuprates has 
proven to a sensitive probe of the pseudogap i"'^^^'^^^'^^'^'^^^ 
Moreover, it is found that for the held perpendicular to 
the layers (which means that Zeeman effects will domi- 
nate orbital magnetoresistance effects) the pseudogap is 
closed at a field given by 



HpG 



(24) 



where T* is the pseudogap temperature. For the hole 
doped cuprates this held is of the order ~ 100 T. In con- 
trast, for the electron-doped cuprates this field is of the 
order ~ 30 T (and T* ~ 30 - 40 K), and so this is much 
more experimentally accessible^ii^ The field and temper- 
ature dependence of the interlayer resistance for several 
superconducting organic charge transfer saltgli^ is qual- 
itatively similar to that for the cuprates. In particular, 
for temperatures less than the zero-field transition tem- 
perature and fields larger than the upper critical field, 
negative magnetoresistance is observed for fields perpen- 
dicular to the layers. A possible explanation is that, as 
in the cuprates, there is a suppression of the density of 
states near the Fermi energy, and the associated pseudo- 
gap decreases with increasing magnetic field. 

A Nernst experiment can be used to probe whether 
there are superconducting fluctuations in the pseudogap 
phase, as has been done in the cupratesJ^ This exper- 
iment is particularly important in understanding the re- 
lation between the pseudogap and superconductivity. 

One could also study the pressure dependence of the 
linear coefficient of heat capacity 7. Since 7 is propor- 
tional to the density of states at the Fermi energy, a de- 
tailed mapping of 7(P) would be an important probe for 
the study the pseudogap. Finally, measurements of the 
Hall effect have also led to important insights into the 
pseudogap of the cuprates^ therefore perhaps the time 
is ripe to revisit these experiments in the organic charge 
transfer salts. 
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APPENDIX A: VERTEX CORRECTIONS AND 
THE DYNAMIC SPIN SUSCEPTIBILITY FOR 
STRONGLY CORRELATED ELECTRONS 

We consider a strongly interacting electron system 
and derive the real and imaginary parts of the dynamic 
susceptibility. We show that under some quite general 
(but specific) conditions that the Korringa ratio is unity. 
Many definitions are simply stated in this appendix since 
most are derived more fully in any number of textbooks 
(for example Ref. 118). The general expression for the 
dynamic susceptibility in Matsubara formalism is given 
by 

Jo 

where /? = l/fc^T is the inverse temperature, r is the 
imaginary time, w„ = (2n + VjirkBT /h is the Matsub- 
ara frequency, ttIq is the component of magnetization in 
the a direction, and Tr is the (imaginary) time ordering 
operator. In order to consider x — h(c[j'J-') we define the 
operators: 



m_(q, r) 



^ P+'i 



cUaM^P.tW' (A2) 



™+(q,T) = ^z^s_^q,TWcpaW- (A3) 

Upon substituting (|A2[) and (|A3[) into (jAl[) and perform- 
ing the appropriate Wick contractions on the operators 
one finds that 

X-+(q, iLOn) = ^ f rfTe^""^r(p + q, -r; p, r) 

xG(p + q,-r)G(p,r)(A4) 

where r(q, t; p, t') is the vertex function, G(p, ipn) is the 
full interacting Green's function given by 



G(p,t) = 



Go(p,r) 



l-Go(p,T)S(p,r)^ 



(A5) 



G"(p,r) is the non interacting Green's function, and 
S(p, t) is the self energy. The r integration can be evalu- 
ated by first transforming the integrand in Eq. (|A4|) into 
momentum space. This gives 



X-+(q, iuJn) 



2/3 ^ 



r(p + q, iprn;P, iPm + «W„) 



X G(p + q, ip„i)G(p, ipm + ii^n) (A6) 



where r(q, ?w„; p, ?w^) is the Fourier transform of 
r(q, -r; p, r') and G(p, ip„) given by 



G(p,zp„) 



iPn - £p - S(p, ip„) ■ 



(A7) 



where Sp is the dispersion of the non-interacting system. 
To evaluate the Matsubara summation, it is convenient 
to express the full interacting Green's function using the 
spectral representation 



G(p,ipn) 



dEi ^(p,^i) 



, 27r ipn - El ' 
where ^^(p, Ei) is the spectral function given by 

-2ImS(p,S) 



(A8) 



{E-ep~ Rel](p, E))^ + (ImS(p, E))^ ' 

(A9) 

Substituting (|A8p into (|A6p . the dynamic susceptibility 
becomes 



X-+(q,za;„) = 2^ 



dEi dE2 



xr(p q, ??3,„; p, ipm + ?t^„) 
^(p-Kq,£;i)A,(p,g2) 

(ipm - Ex)(ipm +iU!n- E2)' 

(AlO) 

At this stage we neglect vertex corrections, that is we 
set r(p -t- q, ipn] p, ipn + iuJn) — 1 for all p, q, p„, and 
LOn- After performing the Matsubara sum and analytical 
continuation iujn w + iri, the dynamic susceptibility is 
given by 



xAs{p,E2 



dEi dE2 

npiEi) - n_F(i?2) 



A,(p + q, El) 

(All) 



where np{E) is the Fermi function. 

First we discuss the imaginary part of x hlli'^)- Us- 
ing the well known relation l/(x-f 177) = P{1 / x) — iTT5(x) , 
where P{y) denotes the principal value, the imaginary 

part of X h(cij*'^n) in the limit of small frequency oj is 

given by 



lim 



(A12) 



The nuclear spin relaxation rate is obtained by summing 
the dynamic susceptibility over all q [c.f., Eq. (Pa|) ] thus 
we find that 



1 



p q 



dnp 
~dE 



(A13) 
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where we have assumed a contact, i.e. momentum in- 
dependent, hyperfine couphng. This expression can be 
written in the more intuitive form 



1 fcs|A|2 dE^^ ( dflF 



where 



(A14) 



(A15) 



is the fuU interacting density of states. At temperatures 
smah enough so that p{E) varies httle with energy within 
ksT of the Fermi energy, Ep, the expression can be fur- 
ther simphfied to 



1 ^ kB\A\^~p^{EF) 



(A16) 



Within the approximation of neglecting vertex correc- 
tions the real part of the dynamic susceptibility is given 
by 

u ^ ^^7? ^ r dE. . 

^ p J-oo 

-^Fip,E,,E2,uj,T) (A17) 



with 



Fip,E,,E2,Lo,T) 



As{p, E2)[nFiEi) ~ nF{E2)] 



huj + El ~ E2 

(A18) 

To perform the integration over E2, we first make the 
change of variable x = lo + Ei — E2 and take the limit 
a; — > 0. Thus the integral over E2 becomes 



J ^F{p,Ei,x,Lo,T) 



—00 
dnp 

JeTi 



By using sum rule 



—As{p, El - x) ^ — ^ - 

Ztt X 



dx 
2^ 



As{p,Ei-x). 



(A19) 



(A20) 



x'_+(q, 0) then becomes 



, , h^Ye r dE^ ^ . .[ dnF\ 

(A21) 

The Knight shift is obtained by setting q = [c.f., Eq. 
icll and is 



\A\le 

27JV 



f°° dE , / dnF\ , . , 



At temperatures sufBciently low that p{Ei) varies lit- 
tle with energy within fc^T near the Fermi energy, the 
Knight shift is given by 



\AWp{Ef) 



(A23) 



Using Eqs. (jA16p and (|A23p . the Korringa ratio for in- 
teracting electrons with a contact hyperfine coupling and 
neglecting vertex corrections is 



/C = 



h 

1. 



IN 

le_ 

IN 



1 



TiTK^^ 

^kB\Alp\EFl 
AttH 



IAMEf) 



(A24) 



For non interacting electrons, the self energy E(q, uj) — 
0; expressions Eqs. (|Al4p and (|A22p are still valid and 
one only need replace p{E) with the non interacting den- 
sity of states p{E) . On the temperature scale over which 
the density of states varies little with energy, the Kor- 
ringa ratio for free electron gas /Cfrco = 1- By comparing 
Eqs. (|A24[) and /Cfrco = 1 we see that any deviation of 
the Korringa ratio from one must be caused by either 
vertex corrections or the q-dependence of the hyperfine 
coupling. 



APPENDIX B: ESTIMATION OF EFFECT OF 
THERMAL EXPANSION OF THE LATTICE ON 
THE KNIGHT SHIFT 

In this appendix we describe how we obtained our esti- 
mates of the isothermal compressibility the linear coeffi- 
cient of thermal expansion and the pressure dependence 
of the Knight shift. By feeding these estimates into Eq. 
(fT8|) we estimate the correction to the Knight shift re- 
quired because measurements are generally performed at 
constant pressure whereas calculations are most natu- 
rally performed at constant volume. The difference be- 
tween these two versions of the Knight shift can be quite 
significant as can be seen in Fig. |4l 

First, let us discuss the first term in the inte- 
grand in Eq. pB]) . We observed that Kf can be 
rewritten as Rp = [h-fl/ {A-KkBl%TiTjqY/'^ . Us- 
ing Mayaffre's data^^ we extracted the pressure de- 
pendence of 1/TiT at constant temperature and esti- 
mated that (l/TiT)^/^ is roughly linear with pressure: 

(l/rir)i/2 3xl0-5p for Ti in second, T in Kelvin, 

and P in bar. We then used this result and the Ko- 
rringa value for non interacting electron gas to cal- 
culate {dKP/dP)j, for K-(ET)2Cu[N(CN)2]Br which is 
found to be around — 3xl0~®/bar. One can compare 
the value obtained here with {dKf /dP)rp obtained from 
the pressure dependence study on effective mass in k- 
(ET)2Cu(NCS)2 compound^^ which is presumably more 
reliable. Since Kg should be proportional to the effective 
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mass, the quantity of interest {dKf /dP)rp can be esti- 
mated from the pressure dependence of the effective mass 
data which yields a value around — 6xl0^^/bar. The two 
estimates agree to with a factor of two. 

Next we need to obtain a value for lattice isother- 
mal compressibility. We use the analytical expression 
obtained from DMFT^ 

(^«r = ^-f^)\e:, (Bl) 

where is the inverse isothermal compressibility 
{vdP/dv)~^, Vq is the reference unit-cell volume, v is 
the unit-cell volume under pressure, Bq is a reference 
bulk elastic modulus, Dq is a reference bandwidth, v is 
a parameter that characterizes the change in the band- 
width under pressure, and Xoi is the electronic suscep- 
tibility We estimated that for k-(ET)2Cu[N(CN)2]C1, 
vo = 1700 1^, Bo = 122 kbar. Do = 0.13 eV, and 
DoXci ^ 1- Putting everything together, the order of 
magnitude of the lattice isothermal compressibility for 
K-(ET)2Cu[N(CN)2]Cl is around -10^ bar. Although 
we are not aware of any measurements of the isothermal 
compressibility of k-(ET)2Cu[N(CN)2]C1 systematic ax- 
ial pressure studiesii^ on a-(BEDT-TTF)NH4Hg(NCS)4 
found that the isothermal compressibility is of order —10^ 
bar, a value which is very close to our crude estimate for 
K-(ET)2Cu[N(CN)2]Cl. 

The temperature dependence of the thermal expan- 
sion at constant pressure has been measured by Miiller 
et al}^^ They found that k-(ET)2Cu[N(CN)2]C1, k- 
(ET)2Cu(NCS)2, and undcuterated and fully deuterated 
K-(ET)2Cu[N(CN)2]Br all have a relatively temperature 
independent thermal expansion above about 80 K while 
complicated features are observed below 80 K associated 
with glassy transitions and the many-body behavior of 
these materials. Since we are only interested in getting 
an order of magnitude, we neglect the complicated tem- 
perature dependence observed in k-(ET)2-X and approx- 
imate the thermal expansion as a constant. The value 
extracted from Miiller et al.'s data is around 10^"* K^^. 



APPENDIX C: ESTIMATION OF ERRORS IN 
FIGURE \6\ 

In this appendix we outline the procedure use to esti- 
mate the errors on the data presented in Fig [B] Let us 



consider a given set of data, for example Tc{P). To a rea- 
sonable degree of accuracy, Tc for k-(ET)2X decreases 
linearly with increasing pressure, i.e. Tc = aP + b where 
a and h are the coefficients obtained from fitting the ex- 
pression to the data. In a typical pressure measurement, 
there will be some uncertainties in the pressure (AP) 
which may be caused by systematic errors due to the un- 
certainties in the pressure calibration. The size of AP 
will depend on a specific apparatus used in the experi- 
ment. For example, a helium gas pressure system would 
have uncertainties as large as 0.1 kbar while a clamped 
pressure cell which uses oil pressure medium may have 
uncertainties as large as 1 kbarii2iiiS2, Knowing AP, we 
can estimate the uncertainty in Tc when it is compared 
with another data set from a different experiment taken 
by a different group, say Tpr^T^{P). This is done by cal- 
culating the following: 



ATc = — ^AP (CI) 



Another way to estimate ATc is to calculate dTc{P)/dP 
from the discontinuity in the thermal expansion and spe- 
cific heat by using the Ehrenfest relatiou;^^ These two 
methods give similar results. For K-(ET)2Cu[N(CN)2]Br, 
dTc{P)/dP is estimated to be around —2.4 K/kbar from 
the first method while it is found to be around —2.2 
K/kbar from the Ehrenfest relation. This procedure is 
repeated for other data sets, Tp^x^{P)i T^y/^{P), and 
7nmr(^') whenever applicable which leads to Tc ± ATc, 
Tp^T^ ±ATp^T^, T/yy/y±AT/\y/y, aud Tnmr i ATnm R ■ 
We then tabulate Pp^T^, ^a-u/d: and Tnmr with respect 
to Tc for a given pressure. The results are shown in Fig 
[51 In some cases we were not able to obtain a reasonable 
fit because either the data are very scattered or there are 
not enough data to perform fit. This is the reason for the 
absence of error bars on the vertical axis for some data 
points in Fig[6l 
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Note that the location of the peak docs not quahtativoly elsewhere then it will simply effect the magnitude of^{Tx). 
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